function ans = p3e4()

	x = linspace(-1,1,200);
	y = f(x);
	
	z = x;
	
	for n = 1:15
		x = [];
		y = [];
		w = [];
		for i = 1:n + 1
			xi = (2*(i - 1)/n) - 1;
			x = [x,xi];
			y = [y,f(xi)];
		endfor
		w = newton(x,y,z);
		hold on
		plot(z,w,'b');
	endfor
	plot(x,y,'r'); 

endfunction

function y = f(x)
	y = 1./(1 + 25*x.^2);
endfunction
